Effective operators from exact many-body renormalization 
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approximation. We analyze the convergence of proton and neutron valence one-body contributions 
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I. INTRODUCTION 

In recent years, ah initio many-body nuclear structure 
calculations, such as the No Core SheU Model (NCSM) 
and Green's Function Monte Carlo (GFMC) have signif- 
icantly progressed, to realistically describe heavier and 
heavier nuclei [J, S, 0, S S I'^ ^he last few years, 
these calculations have been able to reproduce observ- 
ables of light atomic nuclei up to A=14. To deal with 
heavier nuclei [A > 15), even at the current level of 
accessible computing power, it is unavoidable to adopt 
model space restrictions which have to be accompanied 
by proper renormalization of bare NN and NNN interac- 
tions. Significant efforts have been devoted to developing 
the coupled cluster theory with single and double exci- 
tations (CCSD) Q, the importance truncation scheme 
and to recast the ab initio NCSM approach by intro- 
ducing a core and few-body valence clusters [To|. Those 
studies are usually focused on binding energies and nu- 
clear excitation spectra. The electromagnetic and semi- 
leptonic operators, on the other hand, have been studied 
less frequently, and, less is known about their renormal- 
ization. One of the recent studies 3, for example, shows 
that the long-range quadrupole operator undergoes in- 
sufhciently weak renormalization in the two-body cluster 
approximation. This is in contrast to short-range strong 
interactions and short-range operators, which are well- 
renormalized, even in the two-body cluster approxima- 
tion. 

To explore the role of higher-body correlations for 
proper renormalization of the long-range E2 operator, we 
can use the valence cluster expansion (VCE) considered 
in Ref . . Since the effective p-shell interactions con- 
structed in this approximation account exactly for six- 



body correlations, it is also possible to construct the ef- 
fective two-body E2 operator which accounts for those 
six-body cluster correlations. 

In this paper, we present a detailed study of the prop- 
erties of the effective E2 operator in the NCSM formal- 
ism when projected onto a single major shell. The con- 
struction of our effective E2 operator, which acts in a 
Ohft valence space, is achieved as follows. We first per- 
formed a Nj^i^^hil NCSM calculation, for both ^Li and 
^He, using the non-local CD-Bonn potential JJJ. This 
NCSM calculation uses as input the effective interaction 
for ^Li, obtained in the two-body cluster approximation. 
Here, N^^x corresponds to the total oscillator quanta (N) 
above the minimum configuration and varies from 2 to 
16. After a renormalization to the iVmax = space, the 
resulting quadrupole moments and E2 matrix elements 
form the one-body part of the effective E2 operator for 
the p-shell. A N^i^^hn ""Li NCSM calculation is then 
performed for the same range of A^max, as before. After a 
similar renormalization to the A^max = space, the ma- 
trix elements of the E2 operator in this case contain both 
one- and two-body parts. By using the results from the 
^Li and ^He calculations, we are able to subtract the one- 
body contribution from the ^Li E2 matrix elements and 
are left with the pure two-body contribution. This step is 
necessary, as the effective operator will contain two-body 
contributions, even though the bare operator does not. 
We are, thus, able to construct an effective E2 operator 
from the one- and two-body contributions as a function 
of A^max- We demonstrate that the two-body effective E2 
operator can conveniently be parametrized in terms of 
j-dependent one-body effective charges. These effective 
charges account exactly for the one-body contributions 
as well as averaged two-body contributions. 
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II. APPROACH 

A. No Core Shell Model formalism 

The starting point of the NCSM approach is the bare, 
exact A-body Hamihonian with the addition (and later 
subtraction) of the Harmonic OsciUator (HO) potential 

where is the one-body HO Hamiltonian 

and Vij (r2, A) is a bare NN interaction , modified by 
the term introducing A- and fi-dependent corrections to 
partly offset the HO potential present in hj: 

^,,(f7,A)^<^-^(rl-r,)^ (3) 

Another offset occurs for the center of mass (CM) part of 
the HO potential and the CM part of the kinetic energy 
operator, both present in Eq.(IT]), with the addition of a 
Lagrange constraint, AcmHcai, and with Acai chosen 
large and positive. Here, Hcm is the HO hamiltonian 
in the CM coordinates. The eigenvalue problem for the 
exact A-body Hamiltonian ([1]) for A > 3 is very compli- 
cated technically, since an extremely large A-body HO 
basis is required to obtain converged results. The two- 
body cluster (2BC) approximation, consisting of solving 
Eq.(IT]) for the a = 2-body subsystem of A particles, is 
commonly used [l[ to construct effective two-body in- 
teractions for solving the A-body problem in restricted 
model spaces. However, the 2BC approach does not ac- 
count for 3- and higher-body correlations in the effective 
interaction. To include the higher-body effects into the 
effective shell model interaction, we have recently devel- 
oped the two- and three-body valence cluster (2BVC 
and 3BVC) approximation, which includes higher-body 
correlations up to 6- and 7-body, respectively. We gen- 
eralize this technique in order to compute effective shell 
model electromagnetic operators, as described below. 



B. Effective electromagnetic operators 

Let us start with the assumption that we have derived 
a set of effective Hamiltonians H",'^"'"' for different values 
of A^max following the prescription given in ml. Accord- 
ing to the previously adopted notation, the first upper 
index (0) indicates that the effective Hamiltonian is con- 
structed for the Ohfl space ( e.g., for the p-shell). The 



second upper index, iVmax, indicates that the given effec- 
tive Hamiltonian, when diagonalized in the p-shell, ex- 
actly reproduces the low-lying NCSM results obtained in 
the iVijiax^jf^ space. The second lower index, ai, refers 
to the order of cluster expansion employed. Finally, the 
first lower index. A, specifies the mass number of the nu- 
cleus, for which this effective interaction is constructed. 
For simplicity we will omit all these indices below and 
will label the effective Hamiltonian, Tij, corresponding 
to a subset of states with total spin, J. The eigenvec- 
tors of the effective Hamiltonian Tij form the unitary 
transformation Uj, which reduces Ti.j in the Ohfl space 
to diagonal form: 

Ej=UjHjUl. (4) 

This same eigenstate matrix Uj can also be used to cal- 
culate the matrix elements of other effective operators, 
^^^^{Xk; J J'), between eigenstates with spins J and J' 
in the Ohfl space: 

Mta, (Afc; J J') = UjOf^,^ (Afc; JJ')ul, , (5) 

where k is the tensor rank of the operator OA.ai (Afc; J J'); 
and A = -E or M denotes electric or magnetic multipole 
radiation. The required operator mapping procedure 
imposed on the matrix elements of the effective opera- 
tor M'^ai ' calculated with the eigenvectors in the Qhil 
space, is that they are identical to the matrix elements 
^baro q£ -j-j-^g bare one-body operator, obtained from the 
NCSM calculation in the large iVmax^f^ space, i.e., 

Mfa.iXk- J J') = Pa,M\j{Xk- Jj')Pa„ (6) 

where Pa^ is a projector into the ai-body OhQ space. 
A standard NCSM procedure exists for calculating ef- 
fective operators, starting from the bare operator in a 
large space [3]. That is, the matrix M^°['^ is the origi- 
nal bare operator O'^Y^ transformed to the eigenbasis in 
the full Nraax spacc analogous to the one expressed in 
Eq.(5). This technique has been tested in a two-body 
cluster approximation; however, it becomes cumbersome 
in the case of higher-cluster approximations. To over- 
come this problem, we first calculate the bare matrix ele- 
ments, M^^^iXk] J J'), using eigenstates obtained in the 
-^max space. Then the matrix elements of effective opera- 
tors can be determined using the inverse transformation 
to the one given by Eq.®, where the M°^^^^{Xk; J J') 
matrix is replaced by the PaiM'°j^^i{Xk\ JJ')Pai matrix: 

Ot,^{Xk- J J') = UlPa,M)tA{^k- JJ')Pa,Uj', (7) 

according to the definition {i.e., Eq.®) of an effective 
operator. The effective 0'^,^^{Xk;JJ') operator can be 
represented in the standard shell model (SSM) format, 
using the valence cluster expansion (VCE) (similar to 
the VCE for the effective Hamiltonian in [lO|), 

^A^ai = ^AAe + Oa,A^ + 1 + E '^A,A^+ni (8) 

n=2 
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where the upper index, n, stands for n-body part of the 
effective operator in the Ov-body valence cluster (ai — 
Ac+a^)', the first lower index A, for the mass dependence; 
the second lower index Ac + n, for the number of particles 
contributing to corresponding n-body part. The indices 
J J' and Xk are omitted on both sides of the equation for 
simplicity. 



III. EFFECTIVE E2 OPERATOR FOR THE 
P-SHELL 

As an example we consider the E2 operator (i.e., A = 
E and k — 2) for nuclei in the p-shell with an s-shell 
core. Since the spin of the core is zero, there is no core 
contribution for k = 2. So the first term in the VCE, 
given by Eq.([8]), vanishes. Taking ^ = 6, we find, that 
the two-body VCE approximation {i.e., when ai = A = 
6, Ov = 2) is exact: 



Ole{E2), 



(9) 



where the one-body part, i.e., Oq^{E2), is determined 
in terms of the proton and neutron one-body matrix ele- 
ments for the p-shell from the NCSM calculations for ^Li 
and ^He, respectively. To emphasize that Og ^{E2) is an 
one-body operator in the p-shcU space, we can represent 
it in the second-quantization formalism, i.e.. 



Ol,{E2)=J2{^\Oi5im\j)alaj, 



(10) 



where the summation runs over all single-particle states 
considered {i.e., pi/2 andp3/2), {'i\Ol^5{E2)\j) is the ma- 
trix element for these states and o| (oj) is a single- 
particle creation (annihilation) operator. Then the two- 
body part of the effective operator is calculated by 



Ol,{E2)^Of,{E2)-Ol,{E2). 



(11) 



Again, to emphasize its two-body nature, we represent 
it in the second-quantization form: 



Ol,{E2) = lj2i'j\Ol%{E2)\ 



sr)ala'jaras 



0UE2), 



where the one-body operator, Ol^{E2) 



(12) 

is given by 

Eq. pO)) . Utihzing the approach outlined in [l^], we have 
calculated the effective p-shell Hamiltonians for ^Li us- 
ing the 6-body Hamiltonians with A^max = 2, 4, 14 and 
hil. = 20 MeV constructed from the CD-Bonn potential 
[ll|. To perform NCSM calculations we have used the 
specialized version of the shell-model code ANTOINE 
[il[ll], adapted for the NCSM JJ]. The corresponding 
excitation energies of p-shell dominated states and the 
binding energy of ^Li are shown in Fig[Tl as a function 

of A^niax- 

Using the obtained wave functions for different val- 
ues of A^inax, we have calculated the E2 matrix elements 
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FIG. 1: NCSM results for the excitation energies of the 
J"" states and the ground state energy of ®Li, calculated in 
Nmaxti^ spaces with the CD-Bonn potential and h^l — 20 
MeV. The experimental spectra and the ground-state energy 
are shown for comparison. 



M^yF{E2). We have used the following relation. 



between the oscillator length parameter b (measured in 
fm) and harmonic oscillator frequency HQ (measured in 
MeV). Then, using Eq.(I7]), we have determined the ma- 
trix elements of the effective operator Oq%{E2), which 
are listed in the last column of Table [T] for A^max — 14. 

To find the one-body part, 5(^:2), of the E2 effec- 
tive operator we have performed similar NCSM calcula- 
tions for ^Li and '"^He (using the same interaction as for 
the ^Li calculation) and have calculated one-body ma- 
trix elements of the effective E2 operator for protons and 
neutrons, respectively. These are listed in Table [IT] for 
bare (A^max=0) and effective (A^max=12,14 and 16) E2 
operator. Entries in italics represent extrapolated values 
as explained in the text. 

Using these one-body matrix elements of the effective 
one- body E2 operator, 5(£'2), we can determine the 
many-body matrix elements of that operator. Thus, in 
the considered case of two valence nucleons, one can cal- 
culate two-body reduced matrix elements (TBMEs) of 
the one-body E2 operator using the following formula 
[il: 

OTjb; J||Q|.5(^2)||j-jg;J0 _ 

v/(2J+l)(2J' + l) 



{j:\\oUE2)\\j:)St,dfabj'+ (13) 




{jn\Ol5iE2)\\jd)Sa,cfcdj, 
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TABLE I: The values of reduced TBMEs of the one-body 
bare, one- and two-body effective E2 operators for ®Li are 
shown in columns 7,8 and 9, respectively. The values of the 
total reduced TBMEs, Dt are given in the last column. These 
results correspond to A^'max = 14. 
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0.307 


2.662 


3 


1 


1 


1 


2 





-2.070 


-3.331 


-0.537 


-3.868 


3 


1 


3 


3 


2 





0.000 


-0.309 


-0.256 


-0.565 


3 


3 


1 


1 


2 





0.000 


0.000 


-0.072 


-0.072 


3 


3 


3 


3 


2 





-1.463 


-2.254 


-0.590 


-2.844 


1 


1 


1 


1 


1 


1 


0.000 


0.000 


-0.165 


-0.165 


1 


1 
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1 
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3 
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0.000 


0.000 
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3 


1 
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0.000 
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0.000 
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-0.061 
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0.101 
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3 
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-0.869 


1 


3 
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0.000 


0.000 


-0.174 


-0.174 


1 
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3 


2 
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2.535 


4.080 
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4.452 
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3 
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0.000 
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-0.148 


-0.148 


3 


1 
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2 
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1.792 


2.449 
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3.210 
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0.000 
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0.000 
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0.000 
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3 
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3 


3 
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-0.462 
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0.000 
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0.000 


0.000 


-0.039 


-0.039 
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2.738 


4.407 


0.407 


4.814 


3 


1 


3 


1 


2 


2 


-2.738 


-3.741 


-0.731 


-4.473 


3 


1 


3 


3 


2 


2 


0.000 


-0.578 


-0.621 


-1.199 


3 


3 


3 


3 


2 


2 


0.000 


0.000 


0.177 


0.177 


3 


3 


1 


3 


3 


2 


2.449 


3.942 


0.424 


4.366 


3 


3 


3 


1 


3 


2 


0.000 


0.517 


0.425 


0.943 


3 


3 


3 


3 


3 


2 


2.449 


2.921 


0.185 


3.107 


3 


3 


3 


3 


3 


3 


-2.682 


-4.131 


-0.722 


-4.853 



where fabj = (— and the one-body matrix el- 
ements, {jP\\Olr,{E2)\\jP), are listed in Table |TT1 Using 
Eq. (fTS]) , we have calculated the TBMEs of the one-body 
bare (A^max = 0) and effective (for TVmax — 14) E2 oper- 
ators 

D, = {j:3^,;J\\OUE2)\\f,r,;J') 

and have listed them in columns 7 and 8 of Table [H re- 
spectively. Taking into account Eq.®, we calculate the 
TBMEs of the effective two-body E2 operator 

D2 = {j:ji,j\\oi,iE2)\\j:f,-j'), 

which are listed in column 9 of Table |T1 Note, that by 
definition, Di and D2 add to yield the total TBMEs, 
A = J\\0-A{E2)\\j^r,; J'), which we have deter- 

mined using Eq.©. 

Analyzing the results shown in Table IH we identify 
three types of the TBMEs: 

1. The TBMEs which have non-zero values for the 
bare one-body part, i.e., il'i(bare)7^ 0: because 
in the case of the bare operator we have only the 
proton part contributing, this means that those 
TBMEs are allowed according to the one-body se- 
lection rules. The corresponding TBMEs of the 
effective one-body part contain proton as well as 
neutron contributions, if the latter is not forbidden 
by the one-body selection rule for neutron single- 
particle orbitals. The TBMEs of the two-body part 
of the effective E2 operator, D2, always contributes 
constructively to the total TBMEs and is about 
20%, on average, of the effective one-body part, 
Di, in magnitude. 



2. The TBMEs which have zero values for the bare 
one-body part, i.e., £'i(bare) = 0, but non-zero val- 
ues for the effective one-body part, i.e., £'i(eff)7^ 0: 
in this case we have transitions which are forbidden 
according to the one-body selection rule for the pro- 
ton and, subsequently, only the neutron contributes 
to the effective one-body part. The correspond- 
ing two-body parts are of the same order as the 
effective-one-body part and, as in the previous case, 
contribute constructively to the total TBME. 

3. The TBMEs which have zero values for both bare 
one-body part, i.e., Z3i(bare)= 0, and effective one- 
body part, i.e., Z?i(eff)= 0: this means that such 
transitions are forbidden according to the selection 
rules for both the protons and the neutrons. Thus, 
there is only a non-zero two-body part, which is a 
factor of 2-3 less than the two-body part for the 
allowed transitions and is about an order of magni- 
tude smaller than the one-body part for the allowed 
transitions of type 1. 
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TABLE 11: The one-body reduced matrix elements of the effective one-body proton (p = vr) and neutron (p — v) E2 operators, 
Oq ;^{E2), using the CD-Bonn potential with TiQ, = 20 MeV (b=1.44 fm). The calculated one-body and secondary two-body 
effective charges are also shown. 





{f,\\Ol^{E2)\\fA, (efm^) 

\t/ till 1_J,0\ / \ \0 Q i ^ \ / 




A'"max=0 iVmax = 12 iV„,ax = 14 iVmax = 16 A^max ^ OO 




TXV'KV-RV-RV'K V 


3 3 
1 3 


-2.925 0.000 -3.999 -0.508 -4.093 -0.522 -4.162 -0.533 -4.524 -0.553 
2.925 0.000 4.711 0.618 4.847 0.636 4.933 0.650 5.417 -0.664 






3 3 
1 3 


1.000 0.000 1.367 0.174 1.399 0.179 1.422 0.183 1.54 7 0.189 
1.000 0.000 1.610 0.211 1.656 0.217 1.693 0.222 1.852 0.227 






3 3 
f 3 


0.000 0.000 0.256 0.169 0.291 0.161 0.333 0.176 0.632 0.211 
0.000 0.000 0.182 0.209 0.221 0.244 0.255 0.248 0.567 0.294 




e't{-2ja,2jt) 


3 3 
f 3 


1.000 0.000 1.623 0.343 1.690 0.341 1.755 0.359 2.179 O.4OO 
1.000 0.000 1.792 0.420 1.877 0.461 1.948 0.469 2.419 0.521 



A. Effective quadrupole charges 

To quantify the scale of the renormalization of the E2 
operator, it is convenient to use the traditional concept 
of effective quadrupole charges. The effective charges 
are defined as rescaling parameters, which indicate how 
strongly the bare E2 operator has to be enhanced in order 
to reproduce E2 matrix elements calculated in the large 
Nmax^^ model space. 

Thus, it is helpful to define the one-body effective 
charges for specified values of A'^max: 



{jP\\Ol%{E2)U) 
{rJ\Otr{E2)\\j^)-' 



(14) 



where p = t: or v denotes a proton or neutron, respec- 
tively. As one may note from Table HIl and Fig. (2), these 
one-body effective charges are j-dependent and have dif- 
ferent values for protons and neutrons. There is stronger 
renormalization for protons than for neutrons, as mea- 
sured in the magnitude of the shift from their bare values 
(1 and 0, respectively). We also note that the renormal- 
ization for the nondiagonal matrix element is somewhat 
stronger than for the diagonal one. 

To estimate the converged values of the effective 
charges, ecff(iVmax — 00), we fit an exponential plus 
constant (see, e.g., Ref. [ig) to each set of results for 
individual charges as a function of A^max, 



ecff(^max) = aexp (-cA,„ax) + eoff(A„ 



00). (15) 



The uncertainty in the extrapolating functional, assumed 
here as having an exponential form, should not produce 
errors larger than 10 %. We have included the results for 
-^max values from 2 to 14 in the fit. The extrapolated 
converged values are shown in the last two columns of 
Table |TT1 From the Table |TT1 we observe that the average 
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FIG. 2: The one-body effective charges, e'^{2ja,2ji,), for p- 
shell space are shown as a function of increasing model space 

ax . 



magnitude of the proton one-body effective charges for 
^max = 16 appears similar to the standard phenomeno- 
logical value of 1.5. However, the average neutron ef- 
fective charge is around 0.2 which is somewhat smaller 
than the usual phcnomenological value of 0.5. The es- 
timated values of converged proton effective charges are 
only about 10% larger than the corresponding N^^x — 16 
values, while the estimated neutron converged effective 
charges are almost identical to the ones for A^ax = 16, 
which were not included in the exponential fit. 

For the TBMEs of the effective E2 operator of type 1 
and 2, one may model the two body part, D2, in terms 
of one-body matrix elements using relations similar to 
Eq. . In the most general case we would have an addi- 
tional one-body part depending not only on the quantum 
numbers of the single-particle orbitals involved, but also 
on the spin values of initial and final states, J and J'. 
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Alternately, we may introduce secondary (J and J' in- 
dependent) effective charges e2(2ja,2jb), employing the 
following reduction formula for the two-body part: 



^(2J+1)(2J' + 1) 
ej(2j,,2j-,)| ^2 I 0:||Or<=(£;2)|b:)JM/aW' + 

e^(2j.,2j,) I I {jl\\Or%E2)\\jl)5,J,,j. 

We determine the optimal values of the secondary effec- 
tive charges, e2(2ja, 2jf,), by performing a minimiza- 
tion procedure for the following quantity 

ja,b,c,d,J,J' 

where the values of D2 are shown in Table H] and the val- 
ues of D2 are given by Eg. p^ . The optimal values of 
62(2^0, 2jb) are shown in Table HIl As noted above, by us- 
ing Eq. (fT4|) we have calculated the j-dependent one-body 
effective charges, e^(2ja, 2j;,), as a function of increasing 
model-space size, as shown in Fig. [21 

Similar to the case of the one-body E2 matrix elements, 
we employ the relation (|15p to estimate the converged 
values of the secondary effective charges. We take into 
account the result for iVmax values in a range from 2 to 
12. The extrapolated values for iVmax = 16 (in italics) 
are shown in Table [H The estimated converged values 
are shown in the last 2 columns. 

Note that the neutron secondary effective charges 
renormalize weakly and tend to an estimated converged 



value of e^(3,3) = 0.211 and e^(3, 1) = 0.294. The pro- 
ton secondary effective charges, on the other hand, only 
slowly converge to the values of ej(3, 3) = 0.632 and 
62 (3, 1) — 0.567 and, even in the A^max = 16, case consti- 
tute only about 50% of corresponding converged values. 

The total effective charges are given by (2jo, 2jb) — 
e{(2ja,2jb) + e^(2ja,2jb). In this equation, e1{2ja,2jb) 
represents the total effective charge and 6^2-) (2jai 2jfc) 

refers to the one- (two-) body component of the total ef- 
fective charge, respectively. These results are plotted in 
Fig. [3l We see that the two-body component contributes 
a small amount to the one-body component, although it 
is larger for the proton than for the neutron. The esti- 
mated converged values are shown in the last 2 columns 
in the Table H 

The secondary charges enhance the proton one-body 
effective charges by about 25 % and the neutron one-body 
effective charges by about 100 %. We see, therefore, that 
the neutron total effective charges are strongly renormal- 
ized by the secondary effective charges, when compared 
to the proton effective charges. 

Finally, it should be noted that we have discussed the 
averaged calculated effective charges. However, partic- 
ular observables may be sensitive to matrix elements, 
which have been averaged out when calculating sec- 
ondary effective charges. In the next subsection we will 
examine several observables with respect to the two-body 
degrees of freedom of the E2 operator. 

B. Role of two-body components for E2 transitions 

in ''Li 

Here we analyze the nature of the E2 transitions in 
^Li. There are only one-body and two-body contribu- 
tions to total E2 matrix elements in the case of two 
valence particles in the p-shell. Thus, we may deter- 




2*1, T=0 -» 1*1, T=0 



Exp 



CM 3 



8 10 12 14 



FIG. 3: The total effective charges e1{2ja,2jb) = 
e1{2ja,2jb) + e2(2ja,2ji,) are shown as a function of Nmnx- 
The dashed lines indicate the corresponding values of the one- 
body charges at A''niax = 14 (see figure [S]). 



FIG. 4: The one- and two-body contributions, as a function 
of increasing model space size A'max, are shown in grey and 
white, respectively, for the isoscalar transition 2f{T = 0) ^ 
lf{T = 0) matrix element. 
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TABLE IIL The calculated and experimental E2 reduced ma- 
trix elements for ^Li. 



J^,T Jf,T 


{M\E2\\Jf), (efm^) 




iVmax Eq.^ Exp. 
12 14 iVmax ^ oo 


2+,0 l+,0 
3+,0 l+,0 


4.25 4.42 5.59 4.5(1.3) 
4.96 5.16 6.80 8.6(3) 



mine exactly the one-body and two-body (without us- 
ing secondary charges) contents of the E2 matrix ele- 
ments. In Fig. [4] we show the one- and two-body contri- 
butions to the matrix element of the isoscalar transition 
2+(r = 0) ^ 1+{T = 0) as a function of iV,nax- When 
the model space becomes large, the two-body contribu- 
tion to this matrix element becomes a significant con- 
tribution. This indicates that isoscalar transitions are 
renormalized strongly and that higher body correlations 
are essential in accurately calculating the E2 operator. 
This is not the goal of this paper to describe the experi- 
mental data, but for this particular transition, we are able 
to reproduce the experimental result, within experimen- 
tal error, provided that the model space is large enough. 
It is also interesting to estimate a converged value of this 
and other E2 matrix elements using an exponential plus 
constant fit as, we did above for the effective charges. 
We compare the results of the fit and available experi- 
mental data in Table IIIII The results presented in Table 
III indicate that the extrapolated theory gives a reason- 
able estimate of the available experimental data, taking 
into account that theoretical NNN forces have not been 
included in this study. 

In Fig. [5] we show the one- and two-body contribu- 
tions to the matrix element of the isovector transition 
3+(r = 0) 1+(T = 1) as a function of TV^ax- Note 
that in this case the two-body contribution is relatively 



3*1, T=0 -> 1*3,1=1 



FIG. 5: The one- and two-body contributions, as a function of 
increasing model space, are shown in grey and white, respec- 
tively for the isovector transition 3j^(r = 0) — » 1'^{T = 1) 
matrix element. 



small and that the total matrix element renormalizes 
very weakly as the model space increases. Therefore, we 
can see an indication that isovector transitions are not 
strongly-dependent on higher-body correlations. Such 
behavior for the isovector transitions has been noted be- 
fore by the authors in [l7j . 



C. Role of two-body components for the 
quadrupole moment of ®Li 

The quadrupole moment of ®Li is notoriously difficult 
to calculate in the shell-model approach. We will now 
present some insight as to why that may be the case. 
In Fig. [5] we show the one- and two-body components of 
the quadrupole moment as a function of increasing model 
space size. One can immediately draw the interesting 
conclusion that the one- and two-body contributions are 
similar in size, yet have opposite signs relative to each 
other. The two contributions, thus, tend to cancel each 
other to give a small resultant quadrupole moment. The 
calculation in the lAhCl space yields a quadrupole mo- 
ment (3i4rio[l+(T = 0)] = -0.02971e fm^. We were also 
able to calculate the quadrupole moment in the larger, 
16hQ space, where the dimension of the model space ex- 
ceeds 805 million. The obtained value, QiehQ[^~^{T — 
0)] = —0.02969 efm^, is very close to the one for the 
14?if2 space. Note that even in the largest model space 
(16 h^) calculation, our calculated quadrupole moment, 
Qwnnll^iT = 0)] = -0.02969e fm^, is about a factor of 
2.5 times too small, when compared to the experimen- 
tally measured value of Qoxp[l+(r = 0)] = -0.0818 cfm^ 
^IS] . However, the NCSM calculation with the CD-Bonn 
potential for h^l = 13 MeV in the lATi^l space results 
in Qiinn[l^{T = 0)] = -0.04939efm2. The NCSM re- 
sults in the same 14hQ space with the same frequency, 
hn = lS MeV, but with the N3L0 interaction [lg| yield 
a slightly different value, -0.06 efm^. Again, our goal 
here is not to reproduce experiment, but to understand 
the physics of how other physical operators, such as the 
E2 operator, are renormalized due to truncation of the 
model space. In the above case, we conclude that the 
small quadrupole moment for ^Li comes from an almost 
total cancellation between the one- and two-body contri- 
butions, arising from the many-body correlations in the 
five- and six-body clusters, respectively. 



D. Applications to the standard-shell model 

We now turn our attention to testing our effective E2 
operator represented in terms of effective charges. In 
order to do this, we perform a SSM calculation of E2 
transitions for ^Li and ^Li and compare the results to 
the NCSM calculations. In the SSM calculations, we use 
the effective E2 operator created from a NCSM GTi^l cal- 
culation specific to the two-body valence cluster approx- 
imation for each nucleus. We will refer to these effective 
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«Li Q(1\,T=0) 



'Dl 



I 



l+2b 



Qexp=-0 . 08178 (164) efm^ 



2 4 6 8 10 12 14 



This is particularly evident in 



'Li for the 



5- 
21 



and 



21 
the 



21 



transitions and the quadrupole moment of 
state. However, 



moment for the JIl 



in the case of the quadrupole 

1^ and 1^ states, for instance, 
two-body contributions modeled using secondary effec- 
tive charges do not correctly approximate the two-body 
parts of the E2 effective operator. This indicates that 
exact two-body E2 matrix elements have to be used in 
order to obtain consistent information about interference 
of one- and two-body contributions. Furthermore, there 
are higher-body correlations, which have been neglected 
{e.g. the 3-body effective interaction and the 3-body ef- 
fective E2 matrix elements). 



FIG. 6; The quadrupole moment of the ground state for Li 
(l+(r = 0)) is shown in terms of one- and two-body con- 
tributions as a function of increasing model space size. The 
one- and two-body contributions and total quadrupole mo- 
ment are depicted as white, grey and black histograms, re- 
spectively. The experimentally measured quadrupole moment 
[la ] is listed on the figure for comparison. 
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FIG. 7: The quadrupole moment (Q3/2) for the J" 

state and the |^ ^ transition matrix element for ^Li 

using the A7 effective interaction and corresponding effective 
charges. The left most column (criss-crossed filling) refers to 
a SSM calculation using only the one-body effective charges. 
The middle column (zig-zag filling) refers to a SSM calculation 
using the total effective charges. The right most column (solid 
black filling) refers to the full NCSM calculation for the same 
transition at A'^max = 6. 



IV. CONCLUSION 

We have constructed effective E2 operators for p-shell, 
which account exactly for up-to 6-body correlations in 
the lATiil space. Our results indicate that 3- and higher- 
body correlations strongly renormalize the E2 operators, 
enhancing the proton part by about 70% and the neutron 
part by about 40% relative to the bare proton part of the 
E2 operator. Using 1-body valence cluster (IBVC) and 
2-body valence cluster (2BVC) approximations, we have 
decomposed the effective E2 operator into one-body and 
two-body parts, where the effective one-body part ac- 
counts for up-to 5-body correlations in the 16hfl space 
and the effective two-body for residual 6-body correla- 
tions in the lATiil space. We have found that the pro- 
ton two-body part of the effective E2 operator consti- 
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interactions for ^^^^Li as the A7(9) interaction, respec- 
tively. 

The results for ^Li and ^Li are shown in Figs. [7] and 
|S1 respectively. We can see in both figures that us- 
ing only the one-body effective charges is not sufficient 
to reproduce the equivalent NCSM matrix element. The 
secondary effective charges add (in general) a small fi- 
nite contribution to the SSM matrix element, thereby ap- 
proximating the NCSM matrix element more accurately. 



FIG. 8: Quadrupole moments and E2 transitions for ^Li us- 
ing the A9 effective interaction and corresponding effective 
charges. The quadrupole moments are indicated by a "QJ", 
referring to the state with spin J. The E2 transitions are 
labeled by values of a spin for initial and final states. The left 
most column (criss-crossed filling) refers to a SSM calculation 
using only the one-body effective charges. The middle column 
(zig-zag filling) refers to a SSM calculation using the total ef- 
fective charges. The right most column (solid black filling) 
refers to the NCSM calculation for the same transition. 
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tutes on average about 17% of the total proton renor- 
malization, while the neutron two-body renormalization 
is about 50% of the total neutron enhancement. Further- 
more, we noted that the renormalization for the one-body 
nondiagonal matrix element, {p^j2\\E2\\pij2) , is about 
20% stronger than for the diagonal one, (P3/2I |-B2| 1^3/2)- 
We have shown that the two-body part of the effective 
E2 operator can be accounted for reasonably well by in- 
troducing secondary effective quadrupole charges. This 
approximation can be used when the E2 matrix elements 
are on the scale of 1.0 efm^. However, for small matrix 
elements (<0.2 efm^) the effective two-body part has to 
be calculated exactly to reproduce the considerable two- 
body effects. The best illustration of this effect is the 
case of the ^Li ground state, when one-body and two- 
body contributions interfere destructively, resulting in a 
nearly vanishing quadrupole moment. 

We have also shown that our effective E2 operators 
can be used to approximately calculate the quadrupole 
moment and E2 transitions for heavier nuclei in the SSM 
formahsm (see for, e.g.. Fig.® and Fig.®). Such calcu- 
lations are useful in predicting what a NCSM calculation 
for the E2 operator would yield, provided it could be car- 
ried out on a sufficiently large computer. Recall that near 
the center of a major shell, such as the p-shell, the num- 
ber of basis states involved in performing a converged 



NCSM calculation would be computationally expensive. 
However, as we have shown, a SSM calculation is easily 
performed and gives a good estimate of what the NCSM 
result for the E2 operator would be. 
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